First study of the three-gluon static potential in Lattice QCD 
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We estimate the potential energy for a system of three static gluons in Lattice QCD. This is 
relevant for the different models of three-body glueballs have been proposed in the literature, either 
for gluons with a constituent mass, or for massless ones. A Wilson loop adequate to the static 
hybrid three-body system is developed. We study different spacial geometries, to compare the 
starfish model with the triangle model, for the three-gluon potential. We also study two different 
colour structures, symmetric and antisymmetric, and compare the respective static potentials. A 
first simulation is performed in a 24'^ x 48 periodic Lattice, with /3 = 6.2 and o ~ 0.072 fm. 
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I. INTRODUCTION 

We explore, in Lattice QCD, the static potential of 
the three-body glueball system composed of three gluons, 
using Wilson loops. The interest in three-body gluon- 
gluon-gluon systems is increasing in anticipation to the 
future experiments BESIII at IHEP in Bcijin, GLUEX 
at JLab and PANDA at GSI in Darmstadt, dedicated to 
study the mass range of charmoniunr, with a focus in its 
plausible hybrid excitations. Even before the glueballs 
are discovered, the study of two-gluon and three-gluon 
glueballs are respectively relevant to the pomeron [l|, Q 
and to the odderon Q. Thus several models of three- 
gluon models have already started to be developed [3, 0, 

The relevance of computing the static potentials in 
Lattice QCD for 3-gluon models is partly motivated by 
the plausible existence of a constituent nrass for the 
gluon. Several evidences of a gluon effective mass of 
600-1000 MeV, much larger than Aqcd, exist from the 
Lattice QCD gluon propagator in Landau gauge, |9l.[l0t. 
from Schwinger-Dyson and Bogoliubov-Valatin solutions 
for the gluon propagator in Landau gauge |ll| . from 
the analogy of confinenrent in QCD to supercondutiv- 
ity 13 , fr om the Lattice QCD breaking of the adjoint 
string Ij], from the Lattice QCD gluonic excitations of 
the fundamental string [1J| from constituent gluon mod- 
els 0, Ha. \Wl i c omp atible with the Lattice QCD glueball 
spectra 13, [la Hi 1201 , and with the Pomeron trajectory 
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for high energy scattering [ij, l^]- Furthermore, even for 
modelling massless gluons, the knowledge of a static po- 
tential would at least provide one of the components of 
the dynamical potential. For instance, the static quark- 
antiquark potential is frequently applied to light quarks. 
The Wilson loop method was devised to extract, from 
pure-gauge QCD, the static potential for constituent 
quarks and to provide a detailed information on the 
confinement in QCD. In what concerns gluon interac- 
tions, the first Lattice studies were performed by Michael 
[13.[2]| and Bali extended them to other SU(3) represen- 
tations [231. Recently Okiharu and colleagues [23, [23l 
studied for the first time another class of exotic hadrons, 
extending the Wilson loop of three-quark baryons to 
tetraquarks and to pentaquarks. Very recently, Bicudo, 
Cardoso and Oliveira continued the Lattice QCD map- 
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FIG. 1: The starfish-like and triangle-like possible geometries 
for the strings in the static three-gluon system. 



ping of the static potentials for exotic hadrons, with the 
study of the hybrid quark-antiquark-gluon static poten- 
tial, [H,!!!- 

In this paper we study the three-gluon potentials in 
Lattice QCD. We address two novel and important ques- 
tions. Noticing that with three gluons two different 
colour singlets can be constructed, symmetric or anti- 
symmetric, we study whether the respective interactions 
are identical or different. This will be further detailed 
in Section II. Moreover, noticing that a gluon may cou- 
ple to one adjoint string, or to a pair of fundamental 
strings, we study whether the potential is amenable to 
a triangle-shaped triplet of fundamental strings or to 
a starfish-shaped triplet of adjoint strings, as depicted 
in Fig. [TJ A similar discussion on the shape of the 
baryonic strings has been addressed in Lattice QCD 

[13, m, s no, In, us, m, m, 35, m] 

In particular, our study of the hybrid system already 
indicated [2a, [29 that it would be interesting to study 
three-body glueballs, relevant for the odderon problem 
Q. Notice that in Lattice QCD, using the adjoint repre- 
sentation of SU(3), Bali [23 found that the adjoint string 
is compatible with the Casimir scaling, where the Casimir 
invariant Ai • Xj produces for the gg interaction a fac- 
tor 9/4 times larger than the qq interaction. With three 
gluons, a triangle formed by three fundamental strings 
might costs less energy than three adjoint strings with 
a starfish-like geometry, depicted in Fig. [TJ The three- 
gluon potential nray be similar to a sum of three mesonic- 
like quark-antiquark interactions, plus a repulsion act- 
ing only when there is superposition of the fundamental 
strings. This question is also related to the superconduc- 
tor (Type-I versus Type-II) model for confinement, where 



flux tubes repel each other in Type-II superconductors, 
while in Type-I superconductors they attract each other 
and tend to fuse in excited vortices [33] . A first evidence 
of QCD string repulsion was indeed found in our very 
recent study of the hybrid potential [2^, [2y| . The under- 
standing of the three-gluon potential in 3+1 dimension 
Lattice QCD will further clarify our understanding of 
confinement. 

In Section II we derive a class of Wilson Loops ade- 
quate to study the static hybrid potential. This paper is 
mainly analytical, and in Section III we discuss theoreti- 
cally the important questions of the best Wilson Loops to 
distinguish the triangle from the starfish string ground- 
states, and of the differences of the symmetric to anti- 
symmetric potentials. In Section IV we present the first 
results of our numerical Monte-Carlo simulations, and 
conclude. 



II. THREE GLUON WILSON LOOP 

We first construct a wavef unction with three gluons. 
This wavefunction will be the starting point of the Wilson 
Loop. Due to confinement, a hadron, system composed 
of quarks, antiquarks or gluons, must be a colour singlet. 

Each gluon is a state of the adjoint, or octet 8, rep- 
resentation of SU{3). With the tensor product of two 
gluons, different representations of SU{3) can be con- 
structed, 



8(8)8 = le8©8el0el0©27 



(1) 



including a singlet 1 and two octets 8. When we cou- 
ple three gluons, we get not just one colour singlet, but 
two colour singlets (plus many other representations), re- 
sulting from coupling this third octet to each of the two 
octets in the right hand side of eq. ^ , 



8(8)808 



8( 



(2) 



To arrive at the wavefunction for the two colour singlets, 
it is sufficient to study the product of two Gell-Mann 
matrices, since it already produces the relevant colour 
singlet and colour octets resulting from eq. ([T]), 



A"A = o*^" + «/q6cA'' + rfafccA"" 



(3) 



and thus the product of three Gell-Mann matrices al- 
ready produces the two possible colour singlets, that we 
single out in the trace of the product of the three Gell- 
Mann matrices. 



tr {X'^X'X^} = 2ifabc + 2dabc , 



(4) 



and thus the two possible color singlet wavefunctions of 
three gluons are. 
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FIG. 2: Wilsons loop for the ggg potential, (a) for the sym- 
metric colour wavefunction and (b) for the antisymmetric 
colour wavefunction 



where the first combination is anti-symmetric and the 
second is symmetric with respect to the exchange of two 
gluons. 

We build the three-gluon Wilson loop operator inspired 
in the three-quark case of the baryon. In the baryon we 
have a colour singlet wavefunction given by 



and the corresponding Wilson loop is 



'3q 



^ijk^i' j'k' 



(6) 



(7) 



Where X, Y and Z are the elementary paths of the three 
quarks, each composed of the product of successive ele- 
mentary links U starting and ending in wavefunctions of 
the form ©. 

In the three-gluon-glueball case we proceed similarly, 
developing adjoint paths X, Y and Z starting either from 
the symmetric, or antisymmetric, colour singlet wave- 
functions ([5]), as illustrated in Fig. [2 Each adjoint path is 
composed of the product of successive adjoint links, cor- 
responding to gluons, composed of matrices of the SU{3) 
adjoint or octet representation, given in terms of the fun- 
damental representation ones by using the formula 

U^ixf' = irr {X''U,{X)X'' [U,{x)]^} . (8) 

Notice that these adjoint links are unitary matrices, as 
expected by a representation of SU{3) , 

^ U^b i/tbc ^ J2 -Tr[U^X''UX'']Tr[X''U^X''U] 



lTr[U^X''U]Tr[U^X''U] 



6 
where we used the Fierz relation. 



(9) 
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A,^.A2; = 2 (^5uS,k - ls,Ai^ (10) 
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FIG. 3: Graphical version of the Fierz relation, showing that 
when two disconnect paths touch each other at the same point 
where a pair of Gell-Mann matrices is summed in their in- 
dices, this is equivalent to connecting the paths in two differ- 
ent ways, both gauge invariant. 



illustrated in Fig. [3l to contract the A** matrices. 

We now explicitly derive the operator for the three- 
gluon Wilson loop. In the limit of arbitrarily large gluon 
masses, a non-relativistic potential V can be derived from 
the large time behaviour of euclidean time propagators. 
Typically, one has a meson operator O and computes the 
Green function. 



\O{t)O{0)\0) 



exp{-Vt} 



(11) 



for large t. Different types of operators allow the defini- 
tion of different potentials. We can construct the three- 
gluon Wilson loop starting from the gluonic operator. 



Ot,ix) 



J a: 



be 
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(12) 



where the second operator O^ is constructed replacing 
fabc by dabc- In eq. p2p the three octets are situated in 
the same point x. Using the Lattice links to comply with 
gauge invariance, the second operator in eq. (J12p can be 
made non-local to separate the three octet operators. 



Oigix,xl,x2,x'S)^fa 



be 



(13) 



Uf,, (x) • • • Uf,, {x + (ri - l)/ii) g"' {x + nfn) 

J ai 

C/^2 (x)--- Ufi-, (x + {r2 - l)/i2) ff'" [x + ri/12) 

J bi 
- c 

C/^3 (a;) • • • t/^3 {x + (rs - l)/i3) g"^ {x + r^fis) 



where we apply the Lattice QCD prescription of linking 
the fields with links, to maintain the gauge invariance 
of our operator. We also assume the sum over repeated 
indices. The non-relativistic potential requires the com- 
putation of the Green functions present in eq. ([TT|) . As- 
suming that the Gluons are static, and that moreover any 
permutation of gluons is left for the future application of 
the present static potential in constituent gluon models, 
the contraction of the gluon field operators provides ad- 
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FIG. 4: Contractions of the three pairs of Gell-Mann matri- 
ces Graphical resulting from one of the three-gluon wavefunc- 
tions. This shows that the three-gluon Wilson loops, present 
in eqs. p7p and (|18|l . are gauge invariant, because they can 
be written as connected paths of Lattice QCD links U. 



joint temporal links, giving rise to the gluon operator, 



< 



fabc fa'b'c' (14) 

C/pi (x)- ■■ Uf_c^ {x + (ri - l)/ii) 

1/4(0, X + Tijii) ■ ■ ■ Ui{t -l,x + rifii) 

Ulit,x + {n-l)fi,)---Ul^{t,x)y^ X 

Utj.2 (x)--- Uf,^ {x + (7-2 - 1)^2) 

1/4(0, X + r2JJ.2) ■ ■ ■ Uiit ~ l,x + r2/t2) 

_ _ 1 66' 

Ul,it,x + {r2-l)fi2)---Ul^{t,x) X 

C/p3 (x)--- Uf,^ {x + (7-3 - l)/i3) 
1/4(0, X + r^fis) ■ ■ ■ Ui{t -l,x + rsfls) 
UUt,x + irs-l)fi3)---Ul,{t,x)Y' . 



We now translate the adjoint links into quark links. 
This is convenient, both to explicitly show that our Wil- 
son loop is SU{3) gauge invariant, and to arrive at a more 
convenient expression for our computer simulations. So 
let us consider the product of two adjoint links, and apply 




FIG. 5: The antisymmetric three-gluon Wilson loop W-jg ex- 
pressed with patlis of quark-like fundamental U links. 



again the Fierz relation to, say, 
5] [/i^'f/a"' = ^iTr[C/itA''C/iA'']Tr[A^;72A=;72 



2^TiiUi^ X''UiU2\^U2^ 



_^iTr[C/itA'^C/i]Tr[C/2A^C/2t] 



J2UiU2 



(15) 



Thus the product of two adjoint Unks is the adjoint of the 
product of two Unks. Iterating this result to the product 
of an arbitrary number of links, we get that all three 
paths present in eq. (|14p verify. 



u^Ao,x)---ulit,x) 



1 



-Tr 



= X , (16) 

where X is the quark path utilized in the Wilson loop 
for static baryon potentials, corresponding to the gluon 
path X In particular, the Wilson loop in eq. (|14p can be 
decomposed in quark paths X, Y and Z , as in Fig. [21 

wig = fabc /a'6'c'Tr{A''XA'''xt} 

Tr|A''rA'''yt|Tr|A"ZA^'zt| , 

Wig = dabc da'b'c'Tr^X-XX^'X^^ (17) 

Tr|A''yA'''y^}Tr{A"ZA'='z^} , (18) 

extending the three-quark Wilson loop of eq. ^ , replac- 
ing the quark fundamental SU{3) path X, by the gluon 




FIG. 6: The symmetric three-gluon Wilson loop W-fg ex- 
pressed with paths of quark-like fundamental U links. Each 
individual loop is labeled hy a, U. 



adjoint SU{3) path X. We also removed the overall 1/8 
factors since the potentials are independent of the norm 
of the Wilson loops. 

We now proceed to completely translate the results of 
eqs. ((T7)) and (|18|1 into fundamental quark paths. We 
express the eqs. ((T7)) and p^ in terms of correlations of 
the quark paths X, Y, Z only. Noticing, 



1 



fabc = -Tr{(A"A''-A''A'')A^} 

4i 

dabc = iTr{(A''A'' + A^A'^)A'=} , 



(19) 



we replace in eqs. (fT7|) and (fTSll the structure functions 
fabc and dabc by traces of Gell-Mann matrices. Then we 
repeatedly apply the Fierz relation (fTO]) . illustrated in 
Fig. El . 

Subtracting and summing the results of the two differ- 
ent contractions of Fig. 21 we get the contribution of the 
respective symmetric and antisymmetric wavcfunctions 
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FIG. 7: Spatial paths, (a) for an equilateral triangle using the 
vertex of a cube, (b) for an isosceles rect triangle , using the 
vertex of a square. 



where we assumed a sum over repeated indices. Following 
a similiar procedure we also get, 



/afccTr [AMAH'S A^C] - (21) 

= 2iTr{A\Tr{B\Tr{C\ - 2iTr{CBA\ 



Using the resuhs (PD|) and (PT]) we finaUy arrive at the ex- 
pression for the Wilson loop for the antisymmetric colour 
arrangement 



to the three-gluon Wilson loops. 



(20) 



-•.Kj^kl {^^^ - ^°'^)rnn ^jtBlkCn 



= 2i {Tr[CBA] - Tr[ABC]) , 



W^g = +4Tr[XY'']Tr[YZ'i]Tr[ZX'i] (22) 

+4Tr[X^Y]Tr[Y''Z]Tr[Z^X] 
-4Tr[XZ^YX''ZY''] - 4Tr[XY^ ZX^YZ^ , 



depicted in Fig. [5] Using the same techniques for the 
operator for the symmetric colour wavefunction. 



2Tr[ABC] + 2Tr[CBA] - \Tr[A]Tr[BC] - \Tr[B]Tr[CA] - \Tr[C]Tr[AB] + \Tr[A]Tr[B]Tr[C] , 

o o o y 



dabcTr[\<'A]Tr[\''B]Tr[X'C] = 

B 
dabcTrlX" A\^ B]Tr[X'C] = 

C 
dabcTr[X'A\^B\''C] 



2Tr[AC]Tr[B] + 2Tr[BC]Tr[A] + -Tr[AB]Tr[C] - -Tr[ABC] - -Tr[CBA] - -Tr[A]Tr[B]Tr[C] , 



(23) 



2Tr[CBA] + -Tr[ABC] - -Tr[A]Tr[BC] - -Tr[B]Tr[CA] - -Tr[C]Tr[AB] + 2Tr[A]Tr[B]Tr[C] , 



and finaly we get , 



W. 



16, 



16, 



3ff 



4Tr[XY'<ZX'^YZ^]+4Tr[X^ZY^XZ^Y] - —Tr[XY^]Tr[X^Y] - —Tr[Y Z^]Tr[Y'^ Z] 

16 32 

Tr[ZX^]Tr[Z^X] + 4Tr[X^Y]Tr[Y^ ZYTrlZ'i X] + 4Tr\Y'' X]Tr[Z''Y]Tr[X^ Z] + — . (24) 

o o 



r 



The results in terms of quark-like Wilson loops, com- 
posed of fundamental links only, are illustrated in Figs. 
El and El 



III. ANALYTICAL DISCUSSION 



The class of Wilson loops WA and W^ formally de- 
rived in Section II still contain degrees of freedom, that 
we may use to increase the signal to noise ratio. In par- 



ticular, the paths linking the fixed positions xi, X2 and 
X3 of the three gluons, remain to be determined. 

Notice that smearing is a standard technique to in- 
crease the signal to noise ratio of the Wilson loop. The 
smearing [38, 39, 40, 411, [i^, |4^ of the spatial links is a 
technique consisting in repeatedly mixing a link to neigh- 
bour staple-like paths. The resulting mixing is unitarized 
back to a SU{i) matrix. The smearing is expected to 
maximize the signal (of the groundstate) to noise ratio 
when the smearing is comparable to the actual width of 
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FIG. 8: We show the difference V^ - V^ of the three-gluon 
potentials of the two operators Wig and W^g as a function 
of the perimeter p of the respective triangle. For the spatial 
geometry of the loops, we utilize the equilateral triangle of 
Fig. [7ja). The different points plotted correspond to different 
smearings in space and to the same smearing of 1 x 1 in time. 



FIG. 9; We show the three-giuon potentials for the two oper- 
ators( W^g and W^g ) as a function of the perimeter p of the 
respective triangle. We utilize both the equilateral triangle 
and the isosceles rect triangle spatial paths of Fig. [TJa) and 
(6). The results are extracted from 141 SU{3) Lattice QCD 
configurations 24"^ x 48, with the smearing of 50 x 0.2 in space 
and of 1 X 1 in time. 



the QCD confining flux tube. Thus we should try differ- 
ent smearings to arrive at the best signal of the ground- 
state, provided by the exponential decay in eq. pT|) . 

Moreover, the Wilson loops WA and W^g defined in 
Section II depend on the position of the point x, initially 
defined in eq. ([12]) ■ Notice however that the actual static 
potential should not depend on this x point. Possibly, as 
long as we keep fixed the points Xi , X2 and X3 , the spatial 
paths connecting these points could also be arbitrarily 
changed, even if they don't meet in a common point x, 
however this remains to be verified. Importantly, we ex- 
pect that the spatial paths closer to the actual position 
of the strings confining the three gluons will maximize 
the signal to noise ratio. 

In Fig. [7] we show two possible different spatial paths 
linking the points xi, X2 and X3. In this paper, for sim- 
plicity, we use only paths parallel to the lattice grid. In 
Fig. [7] (a) we place the three gluons at the vertices of 
an equilateral triangle, constructed with the egdes of a 
cube. Placing the vertex of the cube at, say (0,0,0), 
three points forming the triangle are, (r, 0,0), (0,r, 0) 
and (0,0, r). In Fig[7I^a), the x point is located at the 
simplest possible position for a numerical simulation, at 
the vertex (0, 0, 0) of the cube. In Fig[7K6), the paths are 
quite simple, we place the three gluons at the vertices of 
an isosceles rect triangle, and the x point coincides with 
the X3, thus the spatial geometry is planar. The paths 
in Fig. IZt^a) and (6) are neither placed at the starfish- 
like string position, nor at the position of the triangle- 
like string position. More sophisticated choices of paths 
might lead to better signal to noise ratios, but the paths 
in Fig. [7] are the simplest for a first simulation. 

On the other hand we may explore analytical similar- 
ities or differences between the Wilson loops W^ and 



W^„. The Casimir scaling, dominating the pertubative 
QCD, and, at least, the short distance potentials, can be 
algebraicly computed. 



(Ai + A2 + As)' 



-(Ai' + Aa' + As'^ 



Ai • A2 = = -6 . 

6 

(25) 

and the result is the same both for the symmetric and 
the antisymmetric potentials. Thus the short range part 
of the interactions should be identical. 

Now we also check that in the limit where two gluons 
are superposed, we recover the normal two-gluon opera- 
tor, where the result is proportional to (the proportion- 
ality factor is irrelevant here). 



VF„, 



W,^ Wq{ 



1 



(26) 



where, say, Wqq =Tr{xyt| jg ^ complete one-quark 
Wilson Loop, to decay exponentially with large times. 
Thus when 2:3 — X2 or equivalently when Z = Y , we get, 
for the antisymmetric loop Wj 



3ff' 



< 



2A{WW* ~1) 



and are also identical in the symmetric loop 1^3^, 



W- 



3g 



40 
-{WW*-1) 



(27) 



(28) 



Importantly, since the result only differs in a physically 
irrelevant constant factor, this shows that whenever two 
of the arms of the starfish are superposed, the two poten- 
tials, for the symmetric and for the antisymmetric cases 
are identical. Then, if any difference occurs, it only oc- 
curs when the arms are separated. Thus we should posi- 
tion the gluons at the vertices of an open triangle, say an 
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FIG. 10: Effective mass plots, for the two operators( W^g and 
wig ), for the equilateral triangle geometry with a perimeter 
p = 15%/2 and for 141 Lattice QCD configurations 24^ x 48. 
In this effective mass plot we use 50 x 0.2 smearing steps in 
space and 1x1 smearing step in time. 



equilateral triangle, or an isosceles rcct triangle, to study 
this possible difference. 



IV. NUMERICAL RESULTS 



where p = ri2 + r23 + r^i is the perimeter of our equi- 
lateral triangle and the sum of the three inter-gluon dis- 
tances, and a is the string tension of the fundamental 
quark-antiquark potential. Both the small difference and 
its linear behaviour confirm our analytical study of Sec- 
tion HI. where we showed that for short distances the 



difference Vs, 



Vasym IS vanishing. 



To verify that the static potentials do not depend on 
the arbitrary meeting point x of the spatial paths, we 
compute the static potentials for two different geome- 
tries, depicted in Fig. [71 In case (a) , the point x is 
placed relatively far from the position of any of the three 
gluons. In case (&) the point x coincides with the position 
of one of the gluons. As anticipated in Section III. the 
potentials show little dependence on the point x. This 
is illustrated in Fig. [21 where both geometries produce 
similar results. 

We also study the absolute size of the antisymmetric 
and symmetric static potentials, defined in Figs. \5\ and 
inland in eqs. (j^^ and (|24p . Again, we explore different 
smearings of the spatial paths. It occurs that the abso- 
lute value of the potentials are more smearing-dependent 
that the nearly smearing-independent difference of the 
potentials. In particular we cannot yet establish precisely 
the strength of the coulomb potential. Nevertheless the 
results of our simulations show in Fig. [SI suggest that 
both the antisymmetric and the symmetric static poten- 
tials are close to the triangle-like model. 



Since this is mainly an analytical paper, in Section IV 
we only numerically simulate the simplest paths to com- 
pute, with the spatial sub-paths depicted in Fig. [3 . We 
perform our simulations with 141 configurations gener- 
ated by the Monte Carlo method in a 24^^ x 48 periodic 
Lattice, with /3 = 6.2 and a - 0.072 fm. 

First we check that the sum of all the different quark- 
like Wilson loops vanish in the limit of large euclidian 
time t. This actually happens, and we also numerically 
check that the Wilson loops , as described in Figs. [Sland 
[6l in the limit of large t tend to. 



h — h — h — h 



tS — tfi 



I7 



1 



(29) 



Then, we study the possible difference between the an- 
tisymmetric and symmetric static potentials, defined in 
Figs. [51 and island in eqs. ^ and ([M]). Notice that we 
explore different smearings of the spatial paths, because 
our spatial paths neither coincide with the vortex posi- 
tions of the starfish-like model, nor with the vortex posi- 
tions of the starfish-like model. Nevertheless the result- 
ing difference, although small, shows little dependence on 
the smearing. The results of our simulations for the dif- 
ference between the antisymmetric and symmetric static 
potentials is show in Fig. [8l suggesting a difference 



i^riangle - 2^-— + '^P' 



(31) 



' V 



Vs, 



Vas 



OMap 



(30) 



where p is the perimeter of the triangle. Our results are 
clearly not compatible with the starfish-like model. In 
the equilateral triangle spatial geometry, the starfish-like 
model corresponds to a linear component of the potential 
component 30 % larger than the one of the triangle-like 
model. The starfish-like a linear term is Icr^mm, where 
Imin is the sum of the distances of the gluons to the Fer- 
mat -Torricelli point, minimizing the total distance of the 
adjoint strings in the starfish-like model. 

Finally, to check that the 24'^ x 48 Lattice configura- 
tions are producing good results, we show the effective 
mass plot of a static potential in Fig. [TUl 

To conclude, we show that there are two, and only two, 
symmetric and antisymmetric, three-gluon static poten- 
tial. We derive the two respective Wilson loops and study 
them analytically. We perform numerical tests, verifying 
that our Wilson loop is correct. Notice that the three- 
gluon Wilson loops include products of up to three fun- 
damental Wilson loops, technically difficult to compute. 
We thus leave the systematic numerical exploration of the 
three-gluon Wilson loops for future works. Nevertheless, 
our numerical simulations already indicate that the sym- 
metric potential is slightly larger than the antisymmetric 
one, and that both are compatible with the triangle-like 
model for the three-gluon static potential. 
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